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We develop a hierarchical structure (HS) analysis for quantitative description of statistical states 
of spatially extended systems. Examples discussed here include an experimental reaction-diffusion 
system with Belousov-Zhabotinsky kinetics, the two-dimensional complex Ginzburg-Landau equa- 
tion and the modified FitzHugh-Nagumon equation, which all show complex dynamics of spirals 
and defects. We demonstrate that the spatial-temporal fluctuation fields in the above mentioned 
systems all display the HS similarity property originally proposed for the study of fully developed 
turbulence [Z.-S. She and E. Leveque, Phys. Rev. Lett. 72, 336 (1994)]. The derived values of a HS 
parameter /3 from experimental and numerical data in various physical regimes exhibit consistent 
trends and characterize the degree of turbulence in the systems near the transition, and the degree 
of heterogeneity of multiple disorders far from the transition. It is suggested that the HS analysis 
offers a useful quantitative description for the complex dynamics of two-dimensional spatiotemporal 
patterns. 

PACS numbers: 05.45.-a; 82.40.Ck; 47.54.-fr; 47.27.-i 



I. INTRODUCTION 



Spatiotemporal chaos have been studied extensively in 
recent years with different level of quantitative charac- 
terization 0. It was discovered that for some spatially 
extended systems such as the Rayleigh-Benard (RB) con- 
vection 1^, an ordered state of straight or weakly 
curved rolls breaks down to a spatiotemporally chaotic 
state consisting of elementary spiral structures which ap- 
pear and disappear in an irregular fashion. Spiral pat- 
terns driven far from equilibrium have also been studied 
in chemical reaction-diffusion systems such as Belousov- 
Zhabotinsky (BZ) reaction Q and in biological systems 
such as cardiac muscle tissue While an individual 
spiral is considered as a self-organized topological defect, 
instabilities can lead to so-called defect-mediated turbu- 
lence that is characterized by an exponential decay 
of correlations with length and time. There have been 
very fruitful efforts in the theoretical, experimental and 
numerical studies on the onset phase of the instability 
which results in the breaking down of spiral waves to the 
early transition to the disorder states (for a review on the 
complex Ginzburg-Landau (CGL) system, see Ref. 0)- 
The study of the transition via instabilities is important 
for problems such as the transition from spiral to spiral 
turbulence in heart tissue which plays an essential role in 
cardiac arrhythmia and fibrillation 

Another question often raised is what quantitative 
measures are useful when spatiotemporal patterns evolve. 
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especially after spatiotemporal chaos are developed 0,I3- 
Attempt to identify defects and conduct statistical study 
of defects in defect-mediated turbulence has yielded some 
interesting results. In the discussion on the CGL sys- 
tem, it was argued that the number of defects can 
be a convenient quantity to characterize spatiotemporal 
chaos. However, the relation between the finite-time di- 
mension and the number of defects has not been fully 
established for systems of all sizes and of all duration 
intervals, and thus the method of separation has not 
appeared fully convincing [?[. In addition, when spa- 
tiotemporal pattern become chaotic, the identification of 
the basic structure (e.g. the defect) is not always possi- 
ble, and its statistics become dubious. In general, spirals 
act more like waves than particles, but sometimes they 
may display spatially intermittent features. Note that, 
although the RB convection, the BZ reaction, and the 
CGL system all exhibit spiral wave behaviors and defect 
structures, the instabilities leading to complex spatiotem- 
poral behaviors are of different nature. Moreover, local 
pattern structures may also appear in various forms of 
spiral, stripe, hexagon, or square. Nevertheless, transi- 
tions from ordered to disordered states exist in all these 
systems, and a common feature of these transitions is 
the generation of multiple scale fluctuations. Therefore, 
simple phenomenology independent of the details of the 
particular system is needed to reveal the essential feature 
of spatiotemporal chaos "buried in the wealth of available 
data" S]. The present work propose such a simple phe- 
nomenology but with a quantitative measure. 

There are several other quantitative measures pro- 
posed for characterizing disorder in pattern dynamics. 
Correlation length was directly estimated in experiments 
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[Pi , and was also extracted from correlation functions lioll 
or from the width of so-called structure factors 0, Il2| . 
Hu et al. proposed to use local wave numbers as 

an order parameter to characterize experimental RB pat- 
terns. Gunaratne et al. [Ij] suggested a disorder function 
to describe disordered patterns close to locally striped 
structures fT5l|. Egolf et al. presented a simple al- 
gorithm for real-time quantitative analysis of local struc- 
ture by computing the local wave number in disordered 
striped patterns. Newell and co-workers used a wave- 
vector field determined by the wavelet transform to study 
the behavior of phase-diffusion equations in the presence 
of defects. All the methods mentioned above are effective 
to describe irregular patterns with "visually" distinct lo- 
cal structures such as local striped structure. When spa- 
tiotemporal chaos are overwhelmingly developed, how- 
ever, local striped structures can not be defined and these 
measures can not be applied. This is the case for spi- 
ral turbulence in highly turbulent regimes of BZ experi- 
ments where either exact defect structures or character- 
istic wave numbers are difficult to obtain. The present 
study gives a different approach which may be applicable 
to much more disordered patterns. 

The BZ reaction system is a particularly interesting 
system for the study of the spiral dynamics, where spirals 
can be supported in excitable or oscillatory media [l|, |j| . 
A recent experimental study (18.] has systematically iden- 
tified several regimes of spatial patterns resulted from two 
instabilities, a Doppler instability and a long wavelength 
instability. At various values of the chemical concentra- 
tion in reservoir A of malonic acid [MA]"^ and in reservoir 
B of sulfuric acid [H2S04]^ , the BZ reaction system dis- 
plays a variety of patterns including (see Fig. 1 which is 
reproduced from Fig. 1 of Ref. jl^ ) 

• simply rotating spiral (S); 

• meandering spiral (Af); 

• chemical turbulence due to the Doppler instability 
(D); 

• abnormal chemical turbulence (A); 

• renascent stable spiral (R); 

• convectively unstable spiral (Ci and C2); 

• chemical turbulence due to the long wavelength in- 
stability (T) 

This rich set of spiral patterns is an ideal system for test- 
ing theoretical methods for characterizing spatiotemporal 
chaos. 

We have proposed a hierarchical structure (HS) de- 
scription of spatiotemporal chaos 18,, 19i] based on an 
earlier HS model of She and Leveque for hydrodynamic 
turbulence poL l2ll | . The starting point of the analysis 
is to construct a multi-scale probability density function 
(PDF) description for the fluctuation field of our inter- 
est, and then to develop a phenomenological model for 
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FIG. 1: The phase diagram of the BZ experiment reproduced 
from [lil. The concentrations [MA]"^ and [H2SOaY^ are the 
control parameters. The solid lines indicate the onsets of 
different instabilities. The dashed lines are the extrapolation 
of the solid lines. The heavy line with arrows at [Myl]'* = 
0.4M indicate the case studied here. 



the description of PDFs over the entire range of scales. 
This set of PDFs carries the probability information for 
fiuctuations at different scales (large versus small) and at 
different intensities (tail events versus core events). The 
scale enters as a parameter when the fluctuation field is 
filtered for defining random variables, and the intensity 
is parameterized by the order of the moments calculated 
with the PDFs. Therefore, the dependence of the mo- 
ments on the scale £ and on the order p constitute a set 
of quantitative measures about the properties of the field. 

The HS model gives a special compact description of 
this set of multi-scale and multi-intensity properties. The 
basic idea in the HS model is that many statistically sta- 
tionary multi-scale fluctuation fields generated from non- 
linear interactions between many degrees of freedom have 
a self-organized property expressible in terms of multi- 
fractal scaling whose multifractal dimension function can 
be understood in terms of a similarity relation between 
structures of increasing intensities of successive moment- 
order p ■ This new similarity relation was proposed 
for turbulence as a generalization of the Kolmogorov's 
complete-scale-similarity |2l|. Later, it was discovered 
that the HS similarity is satisfied in a variety of nonlinear 
systems and other complex systems, such as the RB con- 
vection |25| , the Couette- Taylor flow , flows in rapidly 
rotating annulus |3| ? the climate variations 25.] , the vari- 
ation of nucleotide density along DNA sequences [2^ , the 
diffusion- limited aggregates 27] , a fluctuating luminosity 
field of natural images 28] and several others . More 
interestingly, the previous studies show that the derived 
HS parameters are closely related to global property of 
the turbulent system. For example, in Couette- Taylor 
fiow, one HS parameter gives a signature of the breaking 
down of the Taylor vortex It is thus hoped that 
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the HS analysis would shed light on the properties of a 
fluctuation field displaying spatiotemporal chaos. 

In [Hill, the HS analysis is briefly applied to describe 
the spiral turbulence in the BZ experiment and in the nu- 
merical solutions of the two-dimensional CGL equation. 
It was shown that the HS similarity is accurately verified 
for various experimental and numerical two-dimensional 
fluctuation fields which include both the ordered spirals 
and spatiotemporal chaos. In the study of the BZ exper- 
iment, it is discovered that the measured HS parameter 
(3 shows a transition when the system undergoes a tran- 
sition from one kind of spiral turbulence state to another 
(e.g. from the state D to the state T in Fig. Those 
preliminary findings have inspired the present detailed 
study of the nature of (3 while a more complex sequence 
of transitions are analyzed and more numerical models 
are considered. 

Specifically, we will examine more experimental data 
with transitions between complex states. In Fig. ^ at 
[MA]"^ = 0.4, the patterns vary through M, S, Ci, A, 
C2, R, and T. During the transition, the amount of or- 
der (e.g. spiral structures), the degree of heterogeneity 
(with the mixing of ordered and disordered structures) 
and intermittency vary with parameters. We will show 
that (3 characterizes such transitions near and far from 
the threshold in the experimental BZ system. We sup- 
plement this with a further detailed study of the spi- 
ral dynamics in the CGL equation [2] and in a two- 
component model equation called modified FitzHugh- 
Nagumo (mFN) model [sol l3ll| which has a number of 
applications in biological systems. In all cases, we focus 
on the validity of the HS description and the explanation 
of the dynamics in terms of the degree of heterogeneity 
and of intermittency that the value of /? indicates. We 
believe that the relationship between the values of (3 and 
the dynamics of spatiotemporal patterns is valuable. 

The paper is organized as follows. In Section II, we 
present briefly the HS model with special emphasis on 
the method of the HS similarity test (/3-test) and the 
explanation of the meaning of (3. Section III is devoted 
to the analysis of a sequence of complex transitions in 
the experimental BZ system. In Section IV, we carry out 
a detailed study of the transition of the spiral dynamics 
in the CGL equation and the mFN equation. Section V 
offers a summary and some additional discussion. 



II. THE HS DESCRIPTION 

The HS model has been originally proposed by She 
and Leveque to describe inertial-range multi-scale fluc- 
tuations in a turbulent fluid. The key concept in the 
model is a hierarchy of moment ratios. Let ei denotes a 
certain characteristic variable characterizing fluctuations 
of the physical system at the length scale i. Define the 
pth-order moments Sp{£) = {\ei\P). The HS model intro- 



duce a hierarchy of functions: 



(1) 



each of which has the physical dimension of eg and hence 
describes a certain amplitude of the fluctuations. The 
function Fp{£) depends on the probability density func- 
tions (PDFs) of ee. For a typical set of PDFs of a turbu- 
lent field and of a spatiotemporal field as we study in this 
work, the function Fp(£) increases with p. In other words, 
Fp {£) at higher p is associated with higher intensity fluc- 
tuation structures. We refer Fp{£) to as the pth-order HS 
function. 

An intuitive idea in the HS model is that in a self- 
organized dynamical steady state, there is a similar- 
ity in the dependence of Fp{£) on £ for different p's 
(0 < p < 00). In the case of scaling dependence on £, 
the simplest similarity is that all Fp{£) have the same 
scaling. This corresponds to the Kolmogorov 1941 (K41) 
theory of turbulence in which the large and small scale 
statistics are completely self-similar. For example, when 
ee represents the locally averaged dissipation, then the 
K41 theory predicts that Fp{£) ~ Fo{£) = {ei) ~ £° (the 
mean dissipation is equal to the average energy flux that 
is constant independent of the scale). This complete self- 
similarity law has been demonstrated invalid in experi- 
ments and numerical simulations during the past decade, 
and the dissipation fluctuations have anomalous scalings 
in £. This is so-called intermittency effect 33J. 

She and Leveque postulate that the most intense 
structure Fao{£) = hmp^oo Fp{£) plays a special role, and 
all other fluctuation structures of finite p obey a hierar- 
chical similarity law, namely: 



Fp+i[£) 
F^{£) 



— An 



Fp{£) 

Foo{£) 



(2) 



where Q < (3 < 1, Ap are independent of £. The term 
Foo {£) can be eliminated by considering the ratio 



A,, 



Fp{£) 

Fp{£) \Fp^^{£) 



(3) 



Both sides of Eq. ||2J) can be computed from the empir- 
ical PDFs (or histograms) of eg calculated from an ex- 
perimental or numerical fluctuation field. The linearity 
in the log-log plot of Eq. |j3Jl can be a direct test of the 
validity of Eq. ||2Jl. This will be refer to as the HS sim- 
ilarity test, or the /3-test [23, I3- Figure |2 shows the 
result of applying the /9-test to the data of experimental 
BZ reaction (for details, see below). The linear relation 
is obviously verified, and the slope (3 can be accurately 
estimated. 

The meaning of the parameter (3 can be readily ob- 
tained from the definition Eq. Q [21I . It describes a 
kind of discrepancy between i^p+i and Fp. When (3 ap- 
proaches 1, high p (intense) structures and low p (weak) 



4 



- 




- 




o |}=1.00 [HjSOJ 


■=0.367M 




o p=0.75 [HjSOJ 


°=1.000M 




* |i=0.61 [HjSOJ 


'=1.600M 



0.0001 0.0002 
l°9,.[F(/)/F 



FIG. 2: The HS similarity test for three patterns from the 
experimental BZ system: [7725*04]^ = 0.367M, l.OM, and 
1.6M. The other control parameter is fixed at [MA]'* = 
0.4M. A straight line indicates the validity of the HS sim- 
ilarity. The slope fi is estimated by a least square fitting. For 
clarity, the second and third set of data points are displaced 
vertically up by a suitable amount. 



structures are alike. This may be realized at either ex- 
tremely ordered states or completely homogeneous dis- 
ordered states. In the former, the structures of various 
intensities are strongly correlated. Below, we will report 
observations of this /? for ordered spiral states. The lat- 
ter situation corresponds to the hypothetical field pos- 
tulated by the K41 theory with complete self-similarity; 
unfortunately, it has not yet been observed in real phys- 
ical system. Another extreme case is the limit /3 ^ 0. 
This is the case in which F^(€) stands out as the only 
kind of singular structure which is responsible for the 
physical process. A mathematical model displaying the 
dynamics of one-dimensional shocks, the so-called Burg- 
ers equation, exhibits behavior close to this description, 
where only shocks at isolated points are responsible for 
the energy dissipation and all statistical moments are de- 
termined by the discontinuities across the shocks. 

The study of hydrodynamic turbulence reveals that in 
many realistic turbulent systems, < /? < 1, namely the 
most intense fluctuations do not completely dominate 
and the fluctuation structures of various intensities are 
not completely alike. The fluctuation structures of vari- 
ous intensities are all responsible for the physical process, 
but they are related by the hierarchical similarity law 
Eq. ||2Jl which we believe is a form of the self-organization 
of the system. 

The results reported below on the experimental BZ 
system and numerical spiral dynamics demonstrate that 
when the hierarchical similarity holds, the parameter 
/3 gives a quantitative description of the degree of or- 
der/disorder and of homogeneity/heterogeneity. For /3 
close to one, the system appears to be orderly homoge- 
neous and orderly self-organized; for moderately smaller 
/3 the system contains a mixture of order/disorder and 



appears to be heterogeneous. The smallest /3 indicates 
the appearance of overwhelmingly disordered and inter- 
mittent state with fully developed spatiotemporal chaos. 
The details of these findings are reported below. 



III. ANALYSIS OF EXPERIMENTAL BZ 
SYSTEM 

The experimental BZ system analyzed here has been 
described in detail in Ref. 18]. An extensive phase 
diagram (see Fig. ^) is established in the two dimen- 
sional plane formed by the concentration \lliSO/^ ver- 
sus [MA]'^, in which various pattern dynamical states are 
identified based (essentially) on visual inspection of ex- 
perimental images or luminosity fields. The HS analysis 
is performed on the total variation field, G{x,y), rather 
than directly on the luminosity (image) field f{x,y): 
which is proportional to the chemical concentration: 



0{i,y) 



' dx 



1/2 



(4) 



Then, we construct a set of multi-scale variables to be 
the coarse-grained local variation field, e^, similar to the 
locally averaged energy dissipation in turbulence, which 
is defined as below: 



1 



G{f{x',y'))dx'dy'. 



(5) 



We examine the statistical properties by calculating the 
moments of ei as i varying for increasing control param- 
eter [H2S04]^. The moments Sp{£) = (]£<;]p) are cal- 
culated by a space-average of 256 x 256 pixels and by a 
time-average over a number of images collected at dif- 
ferent times. Furthermore, the HS functions Fp{£) are 
calculated and the ratio across successive orders are eval- 
uated and plotted, as in Fig. |21 It is readily seen that 
the hierarchical similarity is satisfied and the HS param- 
eter (3 can be obtained (by a least square fitting) with 
a good accuracy. Technically speaking, this completes 
the HS analysis of a fluctuation field at a given control 
parameter. 

Detailed examinations of the values of (3 at various con- 
trol parameters constitute a further step in the analysis. 
The variation of /? for all patterns observed in the exper- 
iment at each parameter value of [H2S04]^ is recorded 
in Fig. 13 The values of (3 seem to exhibit strong fluctua- 
tions, and it seems that there is no simple correspondence 
between the value of (3 and the regime. More careful in- 
spection reveals, however, that (3 describes the property 
of the pattern as expected from the theoretical consider- 
ation above. We attempt to itemize our observations as 
follows: 

(a) For the ordered patterns, the values of f3 measured 
are close to one. In Fig. ^a), we plot a typical experi- 
mental image of the ordered spiral state (M) with (3=1. 
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FIG. 3: The HS similarity parameter /3 as a function of 
the control parameter [H2S04]'^ with fixed [MA]"^ = OAM. 
Symbols indicate the pattern regimes as found experimentally 
(see the text for detailed explanation). 
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FIG. 4: Optical images of typical experimental states ob- 
tained with different [H2S04]^ : (a) an ordered state with 
0.367M, (3 — 1.0; (b) a heterogeneous disordered state (co- 
existence of spiral waves and spatiotemporal chaos) with 
0.833M, /3 — 0.7; (c) an abnormal turbulent state with 
0.867M, P = 0.94; (d) a homogeneous intermittent disor- 
dered state with 1.4M, /3 = 0.62. The parameter [MA]"^ are 
fixed at 0.4M. The region shown is 8.5 x 8.5mm^ out of a 
disk with 20mm^ in diameter. 



Theoretically speaking, /3 = 1 implies that low and high 
intensity structures are alike, which happens with an or- 
dered pattern. Note that some values of measured /3 are 
not close to one for ordered spiral states (e.g. the two 
S states for which (3 « 0.75 in Fig.O. One possible ex- 
planation for the exception of the two S states with low 
f3 is the existence of experimental noise (e.g., dark spots 
in the experimental images) that introduces artificially 
sharp structures in the field. Our studies show that such 
"intrusion" of sharp objects lead to a decrease of /3. 

(b) Near the transition between one regime and an- 
other, we observe a depletion of the value of /3 com- 
pared to that measured for the states far from the tran- 
sition (see Fig. EJ. We plot a typical experimental im- 
age near the transition between Ci and A in Fig. EJb) 
with P = 0.7. Theoretically speaking, a lower value of 
P indicates a higher intermittency of the spatial pat- 
terns. Fig.^Jb) indeed shows a mixture of ordered spirals 
and disordered irregular spiral elements, which displays 
clearly a strong heterogeneity. We conclude that one kind 
of intermittency in spatially extended systems is related 
to heterogeneous composition of patterns of different ori- 
gins. In addition, we conclude that a sharp variation of 
f3 usually indicates a transition of regimes. This finding 
is interesting for defining statistical regimes. We believe 
that this is the most interesting finding of this work. 

(c) In fully developed spatiotemporally chaotic states, 
the spatial patterns develop small-scale intermittent 
structures and yield even smaller values of /3 (< 0.65). In 
Fig. ^d) , we show a experimental image for the chemi- 
cal turbulent state T with f3 = 0.62. It is clear that very 
small characteristic scales are developed in the patterns 
which signify large fluctuations in the total variation field 
(because of the spatial derivatives). This enhanced inter- 
mittency appears to be stronger than noted in item (b). 

(d) There are also other disordered states (e.g. A 
and 6*2 ) for which (3 is close to one. Close inspection 
shows that these disordered states all encompass large- 
scale irregular spiral elements, as shown in Fig^Jc) where 
/3 = 0.94. These large-scale patterns are quite homo- 
geneous and spatially ordered, which may yield a large 
value of (3. Note that the HS analysis conducted in this 
paper is a spatially multiscaling analysis which describes 
only the property of instantaneous spatial patterns. Our 
conclusion is that the states in the regimes A and C2 in- 
volve only spatially ordered patterns with actively tem- 
poral dynamics. They are distinguished from the other 
chaotic states like T and D. 

The above analysis helps to understand the nature of 
the parameter (3 in the HS analysis of spatiotemporal 
chaos, and confirms our previous finding that /3's varia- 
tion reveals a transition between two different spatiotem- 
porally chaotic states. The spiral turbulence state due to 
the Doppler instability has a characteristic scale larger 
than that in the spiral turbulence state due to the long 
wavelength instability. Both states are homogeneously 
chaotic, but the later contains more intermittent small- 
scale structures. The HS analysis is suitable to capture 
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FIG. 5: Snapshots of the field f{x, y) from the CGL equation 
at four different parameter values of C3. (a) stable spiral, 
C3 = 0.6; (b) coexistence of spiral and spatiotemporal chaos, 
C3 — 0.74; (c) coexistence of spiral and spatiotemporal chaos, 
ca = 0.75; (d) fully spatiotemporal chaos, ca = 0.8. 
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FIG. 6: The HS similarity test for the four patterns in Fig.|3 
of the CGL equation. A straight line indicates the validity of 
the HS similarity. The slope j3 is estimated by a least square 
fitting. For clarity, the second, the third and the fourth set of 
data points are displaced vertically up by a suitable amount. 



such difference. 



IV. ANALYSIS OF NUMERICAL SPIRAL 
TURBULENCE 

We now turn to the analysis of numerically simulated 
two-dimensional spatiotemporal dynamics. First, in nu- 
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FIG. 7: The HS similarity parameter /3 as a function of the 
control parameter ca in the CGL equation. Note a clear tran- 
sition from an ordered state to a spatiotemporally chaotic 
state. 




FIG. 8: A sequence of evolving fields f{x, y) from the CGL 
equation at difi^erent times, (a) t = 390; (b) t = 500; (c) 
t = 550; (d) t = 770. 



merical calculations, the interaction mechanism respon- 
sible for the generation of spatial patterns is under good 
control, so that spatial patterns are free of noises and 
perturbations that are inevitable in experiments. For ex- 
ample, optical artifacts in the experimental images pose 
great challenge to the HS analysis (in the calculation of 
the total variation field) and thus compromise the clarity 
of the finding. Secondly, the experimental states ana- 
lyzed above are all asymptotical states. We are also in- 
terested in the characterization of transient states. In 
other words, we are interested in testing whether the 
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tion is the two-dimensional CGL equation: 
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FIG. 9: The HS similarity test for the four patterns in Fig.|51 
of the CGL equation. A straight line indicates the validity of 
the HS similarity. The slope /? is estimated by a least square 
fitting. For clarity, the second, the third and the fourth set of 
data points are displaced vertically up by a suitable amount. 
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FIG. 10: The HS similarity parameter /3 as a function of time 
for a sequence of evolving patterns in the CGL equation. Note 
that j3 acts as an order parameter revealing the transition 
from an ordered state to a spatiotemporally chaotic state. 



HS description can characterize the property of a spa- 
tial pattern which undergoes actively dynamical evolu- 
tion. Numerical simulations provide a detailed sequence 
of patterns under dynamic evolution. Finally, the present 
analysis of numerically simulated field may inspire fur- 
ther studies of a variety of two-dimensional partial dif- 
ferential equations of pattern dynamics, e.g. the Navier- 
Stokes equation for the vortex dynamics. Such empirical 
studies are important to the finding of right order pa- 
rameters of chaotic pattern dynamics. 

We will study two model equations for the pattern dy- 
namics leading to spatiotemporal chaos. The first equa- 
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where A{r, t) is a complex function of time (t) and space 
(r) , the real parameters Ci and C3 are coefficients charac- 
terizing linear and nonlinear dispersion respectively. Al- 
though the CGL equation is a normal form relevant only 
near the threshold of a supercritical Hopf bifurcation, 
it has become a popular model to study spatiotempo- 
ral chaos where oscillations and waves are present jjj, Q ■ 
Indeed, our simulation gives rise to both simple spiral 
waves and chaotic states which encompass many disor- 
dered spatiotemporal oscillations. 

The numerical study of Eq. © is performed using an 
Euler algorithm with no-flux boundary condition on a 
256 X 256 square lattice and a time step of At — 0.03. 
The initial conditions are chosen to be a simple spiral 
wave which is the stable solution when ci = —1.40 and 
C3 = 0.50. We then keep ci fixed and increase C3 to 
explore possibly alternative patterns. We choose to an- 
alyze the fluctuation properties of the real part f{x,y) 
of the complex field A{x,y) = f{x,y) + ig{x,y). The 
details of the computation and the transition leading to 
spatiotemporal chaos can be found in Ref. p^. where 
we have established the fact that the /3-test is satisfied 
in both ordered spiral states and spatiotemporal chaos 
states. In the present work, we report the result of a 
systematic HS analysis of a sequence of patterns with 
varying parameters. 

Four snapshots of the field f{x,y) are shown in Fig. \^ 
which contribute to a so-called far-field breaking-up pro- 
cess. As C3 increases, the convective instability occurs 
and the spiral wave breaks up far from the tip. In a range 
of intermediate values of C3, the spiral wave coexists with 
spatiotemporal chaos surrounding it. Finally, at high 
values of C3, the spiral wave totally breaks up and spa- 
tiotemporal chaos dominate the whole box (C3 > 0.76). 

The HS analysis of the solutions of the CGL equation 
begins with the construction of the total variation field. 
Let G{x, y) be defined as in Eq. Q with the field f{x, y). 
The coarse-grained local variation field, e^, is defined as 
in Eq. I^l. The same procedure of the /3-test as explained 
in the previous section is now applied to the data £1. 
As shown in Fig. |H1 good linearity and hence the HS 
similarity are obtained for all the states of both ordered 
spirals and irregular spiral turbulence, consistent with 
the result reported previously and similar to that in 
Fig. 121 Note again that the values of /3 are close to one for 
ordered spiral states, but significantly below one (< 0.9) 
when spatiotemporal chaos are developed. 

The variation of /? as a function of C3 is shown in Fig.^l 
All ordered spiral states have (5 very close to one, and 
the appearance of the irregular patterns surrounding the 
central spiral leads to a smaller /3. It is interesting to 
note that the final chaotic states at high C3 have a /3 w 
0.9 which is bigger than that of the transitional states 
with a mixture of ordered spirals and chaos. This fact 
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FIG. 11: Snapshots of the field of the fast activator u in the 
mFN equation, (a): simple spiral, a = 0.0525; (b): mean- 
dering spiral, a = 0.07; (c): outbreak of breaking-up spiral, 
a = 0.0725; (d): spatiotemporal chaos, a — 0.075. Other 
parameters are a = 0.84, 6 = 0.07. 



is consistent with the previous finding in the analysis of 
the experimental BZ system that the mixed states near 
the transition are more heterogeneous and tend to have 
smaller /3. 

Another interesting aspect of the numerical CGL dy- 
namics is the possibility to follow closely a spontaneous 
generation of spatial disorder from an ordered state. A 
sequence of dynamically evolving patterns arc shown in 
Fig|Sl where the spiral core breaks up without meander- 
ing and then subsequently develops to a complete spa- 
tiotemporal chaotic pattern. This fast breaking-up pro- 
cess is phenomenologically similar to that of the Dopplcr 
instability 0,1121, but with the spiral center fixed. The 
details of the explanation for this process can be found 
in Ref. Here, we show four typical snapshots of the 
spatial pattern at successive times with fixed control pa- 
rameter ci = —1.40 and C3 = 0.8. 

Figure |51 reports the result of the /3-tcst for the corre- 
sponding four pattern states in Fig. |S| Good HS similar- 
ity property is found for all pattern states, including the 
single spiral (/? w 1), the intermediate state with gener- 
ated defects {(3 < 0.9) and the full developed turbulent 
state {(3 « 0.9). Figure ITUl shows the variation of (3 as 
a function of time during the core breaking-up. The pa- 
rameter f3 begins with a value close to one and undergoes 
a transition when disorder is developed. At i « 550, /3 
reaches the minimum value when disorder is mixed with 
order. Later, (3 increases again because of the homogene- 
ity of the disorder. These facts are consistent with the 
observations reported earlier. Once more, the HS param- 
eter /3 acts as an order parameter to the description of 



the evolution of the CGL system from the ordered spiral 
to the disordered spatiotemporal chaos. 

The second model equation we have studied is a sim- 
ple activator-inhibitor model: the two-dimensional mFN 
model [33, 1^. While the CGL equation describes an 
oscillatory dynamics (with two components representing 
the real and imaginary parts of a complex variable field), 
the mFN model is a two-component model of an excitable 
system in which a fast activator (u) and a slow inhibitor 
(v) are in interaction: 



1 , 

Ut = u[u — 1) 

a 



Au, 




vt = fiu) - V, 



6.75u(u- 1)2, 



(7a) 



(7b) 



< It < i 

\<u<l (7c) 
u> \ 



where, the form of f{u) describes an inhibitor produc- 
tion only above a threshold value of u, and the value of 
cr related to the time scales of both activator (w) and 
inhibitor (v) measures the degree of excitability. A nice 
feature of the mFN model is that the degree of the ex- 
citability of the medium can be continuously modified by 
varying the parameters. It is known ^30.] that fast activa- 
tor dynamics (0 < cr ^ 1) and suitable choice for other 
parameters (a < 1, 6 > 0) lead to excitability. When 
one fixes the parameters a and 6 and gradually increases 
CT, the excited spiral waves will loss its stability, begin to 
meander, then break up at the tip, and spread chaotic 
spiral elements outward. Finally, the system develops a 
state of spatiotemporal chaos. This scenario is illustrated 
in Fig. 111! which contains a sequence of images of gray- 
scale picture of the activator u obtained at various times 
of a dynamic evolution. 

The numerical calculation of Eq. ||7J) reported above 
uses an Euler algorithm with no-flux boundary condition 
on a 256 x 256 square lattice. The data presented in 
Fig. ^] are obtained with fixed parameters a = 0.84 and 
h — 0.07, but with several different values of a. The 
simulation starts with a simple spiral initial condition 
that is the stable solution of the system at smaller values 
of a (e.g. a — 0.50). At cr = 0.725, a transition from the 
simple spiral to spatiotemporal chaos as described earlier 
is fully realized. 

To carry out the HS analysis, the total variation field 
G{x, y) for the mFN model is calculated from the fast 
activator field u{x,y) as in Eq. The coarse-grained 
local variation field, e^, is defined as in Eq. 0. We have 
applied the /3-test for the moments of e^, and the result 
is shown in Fig. ^1 Clearly, the hierarchical similarity is 
verified for all values of a and the values of f3 are derived 
accurately in all cases. The variation of /3 as a function 
of cr is presented in Fig. 1131 The correspondence between 
the degrees of orders/disorder and the values of [3 is quite 
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FIG. 12: The HS similarity test for the four pattern states 
in Fig. Illl of the mFN equation. A straight line indicates the 
validity of the HS similarity. The slope /3 is estimated by 
a least square fitting. For clarity, the second, the third and 
the fourth set of data points are displaced vertically up by a 
suitable amount. 



1.00 - 



0.90 - 



0.85 - 




0.80 - 



0.05 0.06 0.07 0.08 0.09 



FIG. 13: The HS similarity parameter /3 as a function of 
the control parameter a in the mFN equation. Note a clear 
transition from an ordered state to a spatiotemporal chaotic 
state. 



clear, as shown in the study of the CGL equation. As 
the parameter a increases, the pattern gradually loses 
the correlation and /3 decreases. Detailed examination 
shows that the tip of the spiral wave begins to meander 
at a = 0.625, and completely breaks up at cr = 0.725. 
The break-up of the tip of the spiral marks the transition 
from /? « 1 to /3 w 0.8. The later state corresponds to a 
fully developed spatiotemporal chaos. 

The above numerical studies reveal that the CGL and 
mFN model equations have less complicated patterns 
than those in the experimental BZ system. This signifi- 
cantly reduces the complexity of the matter, and renders 
the HS description of the dynamics easy to interpret. It 



is worth noting that the numerical data are free from 
complicated disturbance so that the values of /3 show a 
consistent trend of variation with less fluctuations, com- 
pared to the results of Fig. O We therefore find the 
convincing evidence that the HS parameter f3 is like an 
order parameter that describes the transition from an 
ordered spiral state to a homogeneous disordered spiral 
turbulence. 



V. SUMMARY AND DISCUSSION 

One of the most interesting problems in the study of 
the pattern dynamics is to find a quantitative measure 
that describes the global property of the pattern in both 
ordered states and turbulent states. When spatiotempo- 
ral chaos prevail and the effective number of degrees of 
freedom is large, many familiar quantities in the study 
of dynamical systems and chaos become hard to evalu- 
ate. On the other hands, traditional statistical measures 
such as correlation length also become ineffective when 
the pattern presents intermittent features. A new sound 
phenomenology is needed to capture subtle correlations 
across multiple scales and multiple fluctuation intensities. 

We present here a sound phenomenology, the HS 
model, that gives a concise description of a spatiotem- 
poral field by its multifractal scaling properties in terms 
of the HS parameter (3. In this work, we present detailed 
evidence, based on analyzing experimental BZ system 
and numerical CGL and mFN model, that the HS simi- 
larity holds for all the fields (ordered or disordered) and 
the derived parameter /3 describes the strength of the cor- 
relation across scales and across intensities. For /3 close 
to one, the pattern is strongly correlated or ordered; a 
smaller (3 implies a higher degree of heterogeneity and 
a mixture of order and disorder; the smallest /3 is re- 
lated to the intermittency associated with fully devel- 
oped spatiotemporal chaos. In particular, our analysis 
has identified two different disordered states in the BZ 
experiment: the heterogeneous disordered state with a 
mixture of spirals and irregular spiral elements and the 
disordered and intermittent state of fully developed spa- 
tiotemporal chaos. In those cases, (3 provides a useful 
characterization of the patterns in transition. 

Based on the above results, we speculate that the HS 
similarity is an emergent property of self-organized com- 
plex nonlinear multi-scale systems in general. This prop- 
erty, like the maximum entropy property of the ther- 
mal equilibrium state, may be difficult to be derived di- 
rectly from the first principle (such as the Navier-Stokes 
equation for fluids or the Newtonian mechanical law for 
gases). It is thus important to conduct more empirical 
studies for proving/disproving the above conjecture, or 
for assessing its limit of validity. 

For the quantitative analysis of spatiotemporal chaos, 
we believe that it is important to go beyond the sim- 
ple statistics for defects, wave numbers, or other local 
pattern properties, since often they are hard to be de- 
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termined when spatial patterns become overwhelmingly 
irregular. The HS analysis is a global method that takes 
into account the properties of the whole field, but em- 
phasizes at the same time subtle correlations and simi- 
larities among structures at multiple scales and multiple 
intensities. While the previous work 0, ^| establishes 
the validity of the HS model for describing both exper- 
imental and numerical spiral patterns, the present work 
gives detailed evidence that the HS order parameter (/3) 
characterizes the degree of heterogeneity and intermit- 
tency of spatial patterns and is particularly suitable for 
describing the evolution of patterns from ordered states 
to fully developed spatiotemporally chaotic states. This 
is demonstrated in the BZ, CGL and mFN systems. Our 
approach may be easily extended to describe other com- 
plex pattern-formation systems with different unit "cells" 
than spirals, such as strips, hexagons, squares, triangles, 
etc., and this is one of the advantages of the HS analysis. 

To make a concise description for all these statistical 
properties, a phenomenology is required. The HS anal- 
ysis is one of the plausible phenomenological proposals 



that rest on empirical observations (/3-test). The fact 
that the HS similarity is found to be valid for a wide 
class of nonlinear processes indicates the possible univer- 
sal mechanism/principle behind these processes. It is not 
yet clear how this universality should be expressed (e.g. 
in terms of a variational principle like the maximum en- 
tropy and minimum entropy production) , but we believe 
that it deserves special attention in future studies. 
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